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Riemann's minimal surfaces are a complete, embeddable, one-parameter family of minimal sur- 
faces with translational symmetry along one direction. R's infinite number of planar ends are joined 
together by an array of necks, closely matching the morphology of a bicontinuous, lamellar sys- 
tem with pores connecting alternating layers. We demonstrate explicitly that Riemann's minimal 
surfaces are composed of a nonlinear sum of two oppositely- handed helicoids. This description is 
particularly appropriate for describing smectic liquid crystals containing two screw dislocations. 



I. INTRODUCTION 

The connection between statistical mechanics and the 
calculus of variations is both powerful and profound. Be- 
cause geometry and elasticity intimately comingle in soft 
materials, the resulting equations of equilibrium are often 
described or approximated by purely geometrical consid- 
erations; interfaces, in particular, classically adopt the 
morphology of constant mean curvature surfaces, typi- 
cally minimal surfaces. The classic soap film problem 
boils down to minimizing surface area subject to pre- 
scribed boundary conditions. Smectic liquid crystals in- 
troduce an additional layer of complexity to this classic 
problem. Although made of layers, smectics develop peri- 
odic order normal to the layers, while each layer remains 
fluid-like. 

Isolating the topological excitations proves extremely 
powerful in the study of ordered media. Many key physi- 
cal properties may be described entirely by their topolog- 
ical defects. Indeed, in the two-dimensional XY model, 
it is possible to study the energetics and phase transi- 
tions entirely in terms of topological variables. Likewise, 
topological defects form in smectics and, in some descrip- 
tions [1, 2], are the principal players in the melting of the 
smectic phase into the nematic phase. 

Moreover, the effects of surface curvature and topolog- 
ical defects are inextricably linked. While it is simple to 
see how a single topological defect in a surface or family 
of surfaces induces curvature into its host, the converse 
problem is not nearly as transparent. Similarly, the ad- 
dition of further defects into an already curved surface is 
not merely a superposition, nor are the interactions of a 
finite number of defects linear. However, the topology of 
the surfaces themselves provides a framework for study- 
ing the collective effects of curvature and topological de- 
fects. A deep theorem in differential geometry states that 
all minimal surfaces are locally constructed out of pieces 
of helicoids and catenoids [3] . 

An explicit deconstruction of surfaces into their topo- 
logical constituents provides a unique and visual lan- 
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guage with which to describe the properties of many di- 
verse systems. In the case of smectics, we have found 
it useful to construct smectic textures by viewing the 
layers as the Riemann surface of meromorphic func- 
tion on the plane encoding the two-dimensional arrange- 
ment of three-dimensional line defects [4-6]. The one- 
dimensional periodicity arises from the multiple sheets 
formed by logarithmic branch points. Again, great sim- 
plification occurs in those configurations that can be rep- 
resented in this way: the gradient terms of the free energy 
specify harmonic solutions, which are highly amenable to 
this Riemann surface description. 

Lattices of passages and pores often appear in both 
membranes and bicontinuous systems [7, 8]. In smectics, 
this morphology is created by a pair of oppositely-handed 
screw dislocations [9]. Because Riemann's minimal sur- 
faces are constructed from a pair of oppositely-handed 
helicoids [3], it shall serve as our model for a system of 
pores in a smectic. In order to calculate the stability 
of this smectic phase, the exact locations of the screw 
dislocations must be known. The morphology of the dis- 
location cores will also allow us to study the crossover 
between locally helicoidal morphology and a network of 
pores. 

In the following we will review the free energy func- 
tional of the smectic liquid crystal, introduce our con- 
struction for including topological defects in a phase field 
construction, and discuss a new manifest decomposition 
of Riemann's minimal surfaces into a pair of helicoids. 



II. SMECTIC FREE ENERGY, SMECTIC 
TOPOLOGY, AND MINIMAL SURFACES 

As smectics are a one-dimensional crystal of two- 
dimensional layers of liquid with equilibrium spacing 
d = 2ir/q sm , it is apt to describe them with a complex 
order parameter V( x ) = IVtoMI exp (— iq sm $>(x.)). The 
magnitude of the order parameter, |?/>o(x)|, determines 
whether or not smectic order is present. However, all 
subsequent calculations occur deep within the smectic 
phase, where ipo(x) = ipo is constant away from defect 
cores which we treat as cylindrical punctures. The phase 
of the smectic order parameter defines the smectic layers; 
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this may be seen by considering the mass-density wave 
representation of the smectic, where layers are defined by 
level sets of the phase field $(x) = nd, neZ. 

In addition to the standard membrane curvature en- 
ergy terms due to the fluid nature of the smectic layers, 
smectics are also penalized if the layer spacing deviates 
from the equilibrium value. Therefore, the conventional 
form of the rotationally invariant, nonlinear smectic free 
energy is given by 



a pair of oppositely charged helicoids. Like the decompo- 
sition of Scherk's first surface into defects [4] , this allows 
us to consider the defects as the principle players. The 
particular closed form expression for Ricmann's minimal 
surfaces reveals an explicit parametrization the cores of 
the defects. Enforcing the smectic boundary condition - 
flat, equally spaced layers at infinity - while maintain- 
ing the topological defects necessitates curvature in the 
dislocation cores. 



= B J d 3 xl^ 
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III. RIEMANN'S MINIMAL SURFACES 



where B is the bulk modulus, A is the penetration depth 
for splay, and N = V$/|V$| is the unit normal vector 
for the layers [10, 11]. An additional term due to the 
Gaussian curvature may be neglected due to the Gauss- 
Bonnet theorem. 

In smectics the simplest case of a topological defect, 
the screw dislocation, joins together an infinite number 
of flat layers along a straight line, taking the form of a 
hclicoid, $ = z — tan -1 (y/x) = z — Im [log(a; + iy)] . The 
compression energy diverges along the line {x = 0, y — 
0,zeB] [11, 12]. To prevent this divergent contribution, 
the smectic locally melts along the defect line, saving us 
from a topological conundrum at the cost of introduc- 
ing a new set of boundary conditions into the problem. 
While exact solutions to the Euler-Lagrange equations 
often prove elusive, surfaces approximating the real layer 
structure are easily obtained, thanks to the zoo of mini- 
mal surfaces. Choosing a minimal surface with the right 
topology as our starting point, the locations where the 
compression energy diverges are known and we can iden- 
tify the defects. Using one-parameter families of minimal 
surfaces as variational solutions gives us a handle on the 
energetics of these systems. 

Merely specifying the zeroes and poles in a phase field 
completely defines the topology of the resulting surface. 
For regular arrangements of defects in two-dimensions, 
elliptic functions are the building blocks of the textures. 
The properties of elliptic functions vastly simplify both 
the analytic and numeric calculations of the energy. A 
simple example is afforded by Scherk's first surface [13], 
with the multivalued height function 
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The last equality expresses the height function as the 
argument of a trigonometric function [4]. 

Because Ricmann's minimal surfaces have an infinite 
number of planar ends, their boundary conditions and 
asymptotic behaviour make them extremely appealing to 
model the interaction between defects in a smectic. In the 
following, we will provide a precise decomposition of the 
one-parameter family of Ricmann's minimal surfaces into 



Riemann's minimal surfaces are a complete, embedded 
one-parameter family of minimal surfaces [14, 15]. Ric- 
mann's minimal surfaces are foliated by circles and, thus, 
may be described by a phase field of the following form: 



$= [x-a{z)] 2 +[y-f}{z)] 2 -R\z). 
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A minimal surface, by definition, must satisfy the con- 
dition: H = | V • N = for the unit layer normal N. 



A. Weierstrass elliptic functions 



Following Nitsche [16], the phase field must satisfy: 
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from which it follows that R 2 {z) satisfies 
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with d 2 = a 2 + b 2 . Similarly, a'(z) = aR 2 (z) and /3'(z) = 
bR 2 (z). By considering solutions of the form R 2 (z) = 
A f(z)+B, equation (5) becomes the familiar equation for 

the Weierstrass p-function, [p'(z)] = 4 p 3 (z) — 52 p{z) — 
93 = 4 [p(» - ei] [p(z) - e 2 ] [p{z) - e 3 ] , where e 1 + e 2 + 
e 3 = 0. We obtain the following solution for R 2 (z) : 



R 2 (z) = 
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a 2 + b 2 



(6) 



-b 2 + |c 2 ) and 53 



| C (a 2 + 6 2 + |c 2 ) 



where 52 = 4(a , - 

are the Weierstrass invariants. For reference, the half- 
period values are 



e x = p(L0 1 /2) = -- + Va 2 + b 2 +c 2 

2 

e 2 = p(wi/2 + w 3 /2) = - c 



e3 



p(wa/2) 



Va 2 + b 2 +c 2 , 



(7) 



for c < 0. Without loss of generality, we set a — 0, as this 
merely fixes the free rotations about the z— axis. 



FIG. 1: Riemann's minimal surfaces are a complete, one parameter family of minimal surfaces foliated by circles. Due to their 
morphology, they are models for a bicontinuous phase in which a network of pores connects neighbouring regions. In particular, 
Riemann's minimal surface has an infinite number of planar ends, which makes them particularly appealing as a model for 
smectic liquid crystals. 



B. Jacobi elliptic functions 

It behooves us to work with Jacobi elliptic functions 
for the remainder of this paper. Jacobi elliptic functions 
are doubly periodic with real period 2K(m) and imagi- 
nary period AiK(l — to), where m is the square of the 

elliptic modulus and K(m) — J^ 2 dt/yl — m sin 2 (i) is 
the complete elliptic integral of the first kind. To con- 



vert between the two equivalent representations of elliptic 
functions, Weierstrass and Jacobi, we use the identity: 

p(z) = e 3 + (ei - e 3 )ns 2 (z Vei ~ e 3 ,m), (8) 

where m = (e-2 — e 3 )/(ei — e 3 ) . Therefore, in terms of 
Jacobi functions, 
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where a 2 — e\ — e 3 — 2\Jb 2 — c 2 = b/ yjm(l — to), and we have used the identity dn 2 (z, m) + TOsn 2 (z, to) = 1,. The 
integral defining the center of the circles is then 



2/ cn(z, m)dn(^, to) r . -, 

dz ds (z,m) = z(l — m) : ; E am z, to) , to 
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where E(z, m) = J* dt\J 1 — to sin 2 z is the incomplete elliptic integral of the second kind. Note that the integral 
£[am(z,m),m] is composed of a linear part, zE(m)/K(m), and a periodic part, Z[am(z,m),m], In terms of the 
Jacobi elliptic functions, the phase field is 



$ = x 2 



z(l — to) — £(am(z, m), m) + cs(z, TO)dn(z, to) ^ — "^fds (z, m) = 0, 



(11) 



where z = az. Since az is merely a rescaling of the z- 
axis, we have the freedom to choose the constant to be 
fj = 2K(m)/d to ensure that the periodicity of the phase 
field matches the ideal spacing of the smectic layers d. 
The final parameter, b = AK(m) 2 ^m{\ — to), may now 



be expressed in terms of the elliptic modulus. The re- 
sulting surface is now a one-parameter family of surfaces, 
depending on to € [0, 1]. 

As there are an infinite number of phase fields that 
describe the same surface, determining the appropriate 
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representation of the phase field for each given applica- 
tion requires careful consideration. First, we demonstrate 
that Riemann's minimal surfaces are composed topolog- 
ically of two screw dislocations of opposite handedness, 
and, secondly, we find a form of the phase field for Rie- 
mann's minimal surfaces that satisfies the boundary con- 
ditions of a smectic liquid crystal. 

Expanding the square in y in the phase field, we obtain 

$ = x 2 + \y — ((z, to)] 2 — r] 2 (z, m) 

+ 2^[y-C(z,m)] cs(z,m)dn(z,m) = 0, (12) 



where we define C(z,to) = (a/b)[z(l — m) — 
E(am(z,m),m)], and f] 2 {z,m) = ?2 [ds 2 (z, m) — 
cs 2 (z, m)dn 2 (z, to)] = (cr 2 /6 2 )dn 2 (2;, to). After some al- 
gebra, we find the equivalent phase field, 



$ = 1 + 



[y- C(z,m)] 2 

x 2 — rj 2 (z, to) 



y-((z,m) y-((z,m) 



■ + rj(z,m) x~r]{z,m) 



ocs(z, ?n)dn(z, m) 
brj(z, m) 
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(13) 



This provides a relation between y, z, and x, 

sc(z,to), (14) 



y-Q(z,m) 



y-((z,m) 
x — 'q(z,m) 



x+r](z,m) x—r](z,m) 

where we note that the definition of sc(z, to) = 
sin (am(z, to)) / cos (am(z, to)) . Taking the arctangent of 
both sides and rearranging the terms yields the phase 
field, 



$ = am((jz, to) + tan 



tan 



y-((crz,m) 
x + rj{az, to) 



y - C(vz,m) 
x — rj(az, to) 

= 0, 



(15) 



which is clearly a nonlinear sum of two helicoids of op- 
posite handedness. The cores of the defects are located 



at x = ±r)(az,m), y = (,{az,m), where ((az,m) con- 
sists of a linear component with a periodic term su- 
perimposed upon it, and r)(az,m) is a periodic func- 
tion with maximum at T](0, to) = f and minimum at 
r](K(m),m) = ^\/l — m. 

In order to satisfy the smectic boundary condition, flat 
evenly spaced layers infinitely far from the defects, the 
phase field must be normalized such that the compres- 
sion energy vanishes at infinity, i.e. lim (V$) 2 = 1. 

x.y— 7-oo 

Since the derivative d z am(z 1 m) = dn(z,m) is not in- 
dependent of z, this phase field does not satisfy this 
boundary condition. One final manipulation, based on 
the identity f (am(z, to), to) = z, where F(z,m) = 
Jq dt(l — to sin 2 i) -1 / 2 is the incomplete elliptic integral 
of the first kind, gives the phase field representation of 
Riemann's minimal surfaces, 



<f> = z ■ 



F 



tan 



' y - C{^z,m)" 
^x - r]((TZ,m) i 
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'y- C,{pz,™Y 
x + r\{az,m) i 



|, TO 



= 0. 



(16) 



This is our main result. 



IV. PORES AND DEFECT CORES 

The entire one-parameter family of Riemann's mini- 
mal surfaces are parametrized by the elliptic modulus, 
to e [0,1]. As it sweeps through all allowed values, the 
morphology of the surface, and simultaneously the shape 
of the defects, changes. In the limit that to — > 0, the de- 



fect cores are straight and infinitely far apart. However, 
as to increases, the distance of closest approach between 
the two defect cores decreases. The elliptic modulus is a 
proxy for the distance between the two helicoids, where 
the maximum, minimum and average separation of the 
defects are respectively d max = -4=, d min = — * 

°v m (7-v/m(l— m) 

and d av „ = , = , shown in Fig. 2. If the layers are 

cr 2 ^/m(l — m) 

to remain flat along the boundary at infinity, the cores 
necessarily bend. Conversely, two straight screw disloca- 
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tions in a smectic will force the layers to bend so that 
H 7^ 0. The amplitude of the distortion of the cores also 
must increase as m increases. Nevertheless, the topolog- 
ical character of straight screw dislocations is preserved 
even though the cores of the dislocations possess cur- 
vature. The line integral of the director field around a 
closed loop encircling a defect, S N • dl — ±d, remains 
constant regardless of the position or orientation of the 




FIG. 2: Due to the curvature of the defect cores, the elliptic 
modulus parametrizes the distance between the two helicoids. 
The maximum, minimum, and average separation are given 

by dmax = /— , ^min = ; = , and rf ave = 



respectively. 



T^J m(l — m) 3 



2 ^/ m(l — m) 5 



Alternatively, we might have used the construction <I> 
via foliation of circles implies that layers flatten when 
the radius approaches infinity and that adjoining layers 
are connected by necks. By varying the surface param- 
eter through the range m g [0,1], the family of Rie- 
mann's minimal surfaces transforms between these two 
morphologies, as seen in FIG. 3. 



(a) 




(b) 





V. ENERGETICS 

Comparing the energetic cost of Riemann's minimal 
surfaces to that of other smectic textures determines their 
relative stabilities. To illustrate how we do this, we first 
consider the case of a single screw dislocation. For such 
a comparison to be valid, all of the divergences must be 
treated identically. While it is commonplace to cut off 
the integration in cylindrical regions, where varying the 
radius determines the functional form of the divergence, 
the curved nature of the defects precluded straightfor- 
ward implementation of this method. Instead, we opt 
for a more physically-motivated cutoff: the defect core is 
identified with the volume in which the energy density ex- 
ceeds a critical value eo. This cutoff captures the fact that 
the real smectic melts when the local energy density ex- 
ceeds the condensation energy for forming a smectic. The 
single dislocation, <£> S crcw = z — (l/7r)tan _1 (y/x), has 

compression energy density / scre -» 

v 2 ) 2 



= I7v$ ) 2 - 

[V v ^ screw } 



H 2 = 



r 4 (x 2 



(nr) 4 . The energy density exceeds 



l/(7rey 4 ), corresponding to a free energy 



V 

eo for r < ro 
per unit length of F st 

In order to calculate the free energy of Riemann's 
minimal surfaces, we numerically integrated F — 
J dVf comp 9(e - /comp), where / comp is the compression 
energy density (B/2)(l — V^ 2 ) 2 and 9(x) is the Heaviside 
theta function, defined to be 1 when x > and when 
x < 0. The resulting energies depend on both the elliptic 
modulus m and the free energy cutoff eo- As with the 
single screw dislocation, we expect the energy to scale as 
the square root of the cutoff, and, indeed, for fixed elliptic 
modulus, the energy fits F m (eo) = A + B^/eo, where the 
m-dependence of the parameters is A(m) = aim/ (I — m) 
and B{m) = b\ + b-ifn/ (1 — m). Thus, the free energy for 
Riemann's minimal surfaces is fit by 

F(m, eo) = a— == + ^e^(b + c- , (17) 

yl — m V 1 — mJ 

for parameter values a = 21.7622 ± 0.0300, b = 0.6319 ± 
0.0009 and c = 0.1387±0.0003. Though there is no justi- 
fication for this functional form, the fit, shown in FIG. 4, 
is quite good. It would be interesting to understand the 
success of this fit. Note that the constant proportional to 
y/eg is close to 2/7T w 0.6366, giving the contribution from 
two non-interacting dislocations. The m dependence be- 
comes a 1/r 2 energetic interaction between the defects 
for large separations r. 



VI. DISCUSSION 



FIG. 3: Two of Riemann's minimal surfaces are shown here 
for parameter values (a) m = 0.075 and (b) m = 0.925. 



Nowhere is the "duality" between the two topological 
descriptions of Riemann's minimal surfaces more impor- 
tant than in the energetics. By choosing the helicoidal 
description of the phase field, the divergences in the ener- 
getics lie along the defect cores. Here the smectic layers 
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FIG. 4: The results of the numeric integration of the free 
energy functional for Riemann's minimal surfaces for varying 
values of the elliptic modulus m and the cutoff energy eo are 
plotted along with the fit in Eq. (17). 

are flat and uniformly spaced at infinity. Conversely, a 
smectic containing a pore will have curvature singulari- 
ties down the center of the pore and at infinity, reminis- 



cent of the focal lines in the focal conic texture. 

This simple example highlights the importance of topo- 
logical decomposition of minimal surfaces in studying 
complex systems in nature. Enumerating all topologi- 
cal defects in a surface simultaneously identifies all en- 
ergetic singularities. This technique will help explain a 
variety of phenomena - in particular, triply periodic min- 
imal surfaces which appear in bicontinuous cubic phases 
of systems ranging from mitochondrial membranes to the 
dark conglomerate phase of bent core liquid crystals to 
binary metallic alloys. The network of defects recovered 
from a topological decomposition should explain both 
their extraordinary stability and complex phase diagram. 
By tuning the interfacial preference for negative Gaus- 
sian curvature, the cubic phases continuously transform 
into one another following the progression Schwarz P sur- 
face to diamond surface to gyroid. Yet simple curvature 
considerations cannot explain this series of surfaces. Al- 
though these surfaces are topologically distinct, tracking 
the evolution of the lattice of their defects may eluci- 
date the means by which these surfaces smoothly change 
topology [18]. 
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